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I Introduction 

We have recently completed the NASA funded study, "Mesoscale monitoring of the soil 
freeze/thaw boundary from orbital microwave radiometry." While most of our objectives were 
met, a few will have to be addressed within the newly begun NASA project, "Mapping regional 
freeze/thaw patterns with satellite microwave radiometry." The salient scientific products of the 
initial project are that regional freeze/thaw maps can be extracted from Nimbus 7, Scanning 
Multichannel Microwave Radiometer (SMMR) data, that diurnal thermal gradients have a small but 
measurable effect upon the SMMR spectral gradient, and that scale-space filtering can be used to 
improve the spatial resolution of a freeze/thaw classified image. 

This Final Report of the initial project contains a summary of the investigation, a discussion 
of accomplishments and of unresolved issues, and copies of publications and of conference 
abstracts and proceedings. We very much appreciate the opportunity to have worked upon this 
exciting, and, we believe, fruitful project. 


II Summary of the Investigation 

Soil moisture contributes to the energy exchange between the air and ground through latent 
heats of fusion and vaporization, and to rainfall runoff through the field moisture deficiency of a 
drainage basin. Whether or not soil moisture is frozen affects both the rate of energy transfer to the 
atmosphere, and the infiltration capacity of the soil. The amount and state of soil moisture are 
regional parameters that one would like to estimate using satellite remote sensing. There is a large 
body of literature about estimating soil moisture from radiobrightness (e.g. Burke et al., 1979; 
Wang et al., 1982; Schmugge, 1983; Jackson et al., 1984; Camillo and Schmugge, 1984; and 
Schmugge et al., 1986). We have developed a technique for mapping the spatial extent of frozen 
soils from the spectral characteristics of the 10.7-37 GHz radiobrightness. Through computational 
models for the spectral radiobrightness of diurnally heated, freezing soils, we identified a 
distinctive radiobrightness signature for frozen soils, and cast that signature as a discriminant for 
unsupervised classification. Our initial results were reported at the Interdisciplinary Science Land 
Surface Climatology Program's (ISLSCP) meeting in Las Cruces, NM, November 16-18, 1988 
by A. W. England, and will appear as a paper in a special NASA ISLSCP journal publication 
(Zuemdorfer et al., 1989). 

Freezing influences the apparent radiobrightness temperature of the ground, T 5 , through 
parameters in the approximation (Ulaby et al., 1981), 

Tb — e T 0 + (1-e) T^y, 

where e and T 0 are the emissivity and surface temperature of the ground, respectively, and T s ^y is 
the effective sky brightness. In this approximation, atmospheric transmissivity is ignored. Frozen 
ground exhibits signatures of (1) lower thermal temperature, T 0 , (2) higher emissivity, e, and, (3) 
a decrease in brightness temperature with microwave frequency, /, 


3Tb 

3 / 


< 0.0 


Signatures (1) and (2) are frequently ambiguous indicators of frozen ground because the 
variations in radiobrightness that result from freezing are easily confused with the effects of 
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variations in soil moisture. Water molecules in frozen plants and soils are not free to align 
themselves with microwave electric fields. This constraint upon the rotational freedom of water 
gives rise to an apparent dryness to microwaves. The consequence is a decrease in the real part of 
the dielectric constant, and an increase in frozen soil emissivity. For example, the real part of 

dielectric constants, e', and corresponding emissivities at nadir, e(0), of two, homogeneous, 

smooth surfaced, 15% moist soils at 10 GHz are (e’ from Hoekstra and Delaney, 1974): 


+ 5° C - 5° C 


Material 

£ 

e(0) 

lb 

£ 

effil 

lb 

Goodrich Clay 

8.2 

0.77 

221 

4.9 

0.86 

235 

Fairbanks Silt 

9.6 

0.74 

214 

4.1 

0.89 

242 


Because of increasing emissivity with freezing, a 10° decrease in the clay and silt soil 
temperatures, from +5° C to -5° C, would cause an increase in T^ of approximately +14 K and 
+28 K, respectively. The positive direction of change in T 5 with soil freezing will cause confusion 
in discrimination between moist soils which will appear radiometrically warmer when frozen, and 
dry soils which undergo little molecular change and will appear radiometrically colder. 

The shift in emissivity with freezing is most pronounced at the lower microwave 
frequencies. At 37 GHz, the effect is reduced but not absent. We observe that the 37 GHz 
radiobrightness correlates relatively well with air temperature (Figure 1). Since soil surface 
temperature should follow the air temperature, the 37 GHz radiobrightness can be expected to 
provide a reasonably reliable estimate of soil surface temperature. However, discrimination based 
only on the 37 GHz radiobrightness would misclassify too often. 

The third signature of frozen soil occurs because freezing reduces the imaginary part of the 
dielectric constant, e", proportionally more than it does the real pan, e'. The loss tangent. 


tan 8=e"/e' 


is a measure of the attenuation per microwave wavelength in emitting media. Reduced loss 
tangent, or lower attenuation, means that thermally emitted photons originate deeper within 
emitting media. That is, the effective depth of emission, Zq, (1-e‘l of the emission originates 

above z e ) becomes a larger fraction of the free-space wavelength, X 0 (England, 1974, 1975, 1976, 
and 1977). For example, Goodrich Clay and Fairbanks Silt exhibit an increase of z e with freezing 
(dielectric data from Hoekstra and Delaney ,1974) , 


+ 5°C - 5° C 



£l 

£ 

tan S 

Ze 

£ 

£ 

tan 5 

Ze 

Goodrich Clay 

8.2 

3.5 

0.43 


4.9 

1.0 

0.20 

0.36 X 0 

Fairbanks Silt 

9.6 

5.0 

0.52 

0.10 

4.1 

0.02 

0.005 

15.7 X Q 
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The effective emission depth of moist soils is typically 10% of the free-space wavelength. 
Frozen soils have effective emission depths that may be 30% or more of free-space wavelength. 
The effective emission depth of frozen sandy soils, like the Fairbanks Silt, can be several 
wavelengths. In the more transparent emitting media, particularly in frozen sandy soil or dry snow, 
the greater average thermal photon path lengths have two effects: (1) a greater likelihood that 
thermal gradients affect spectral gradients, and (2) a greater opportunity for volume scattering of 
photons. 

(1) Thermally induced spectral gradients occur because longer wavelength photons tend to 
originate below the optical surface where thermal temperatures may differ by several degrees from 
surface temperatures. For the lower loss tangents of frozen soil, this difference in average emitting 
depth is enough to reflect near surface thermal gradients caused by diurnal heating. That is, a 
positive thermal gradient, 3Tq/3z, will yield a negative spectral gradient, 3T^/3/. SMMR data are 
collected at midnight and noon. In the absence of changing weather conditions, midnight thermal 
gradients will be positive and noon thermal gradients will be negative (Figure 2) so that midnight 
spectral gradients will be negative, and noon spectral gradients will be positive. An average +0.2 
Kelvin//(GHz) shift in the spectral gradient is observed between midnight and noon for SMMR 
radiometric brightnesses (Figure 3). While we have developed a computer model of these gradient 
effects, for the purposes of this report, thermally induced spectral gradients are noise to be filtered 
out. 


(2) The second consequence of soil freezing is a greater opportunity for volume scattering - 
- particularly at shorter microwave wavelengths. This occurs because of the greater average 
photon path lengths in frozen soil, and because plants and soil appear increasingly heterogeneous 
at shorter wavelengths. This "law of darkening" means that, for an isothermal, volume scattering 
halfspace, 


3Tb 

3 / < 


0.0 


(England, 1974). Frozen terrain may also be snow covered. Dry snow is exceedingly transparent 
to microwaves so that snow exhibits significant of darkening (Figure 4, Edgerton et al., 1971). 
That is, both frozen soil and snow tend to exhibit negative spectral gradients. While neither a low 
37 GHz radiobrightness nor a negative spectral gradient is solely adequate as a classifier of frozen 
soils, particularly at the relatively coarse resolutions of the Nimbus-7 SMMR, a discriminant based 
upon a combination of these signatures appears to classify correctly most of the time. 

SMMR radiobrightness data at 6.6 GHz, 10.7 GHz, 18 GHz, and 37 GHz were obtained 
for August 1, 1984, through December 31, 1984, over an area that included North Dakota, about 
half of each neighboring state, and part of southern Canada (Figures 7-9). We chose this large, 
relatively uniform area because of the low spatial resolution of the SMMR instruments — 150 Km 
at 6.6 GHz, 100 Km at 10.7 GHz, 60 Km at 18 GHz, and 30 Km at 37 GHz, and because of the 
importance of soil moisture state to this region’s hydrologic processes. The data arrived from the 
National Space Science Data Center (NSSDC) on 21, high density, SMMR Cell Tapes. Such data 
are referenced to latitude and longitude in a satellite-centered coordinate system. We produced two 
types of image products: Single-band, radiobrightness images at the intrinsic resolution of each 
sensor, and (2) composite, multi-band images at a common resolution based upon local area 
averaging. Each radiobrightness pixel was referenced to latitude-longitude in a Mercator projection 
by interpolation and resampling the Cell Tape data. We used a bi-cubic approximation of a sine 
function (Moik, 1980) for the interpolation. H and V radiobrightnesses were averaged to produce 
a single brightness for each pixel for each frequency. 
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In addition to large area images, local area spatial averages of radiobrightness were 
calculated for each radiobrightness channel at 7 meteorologic sites within our test region-Miles 
City, MT; Bismark, Fargo, and Williston, ND; and Abileen, Huron, and Rapid City, SD. A local 
area is defined as a 150x150 Km cell centered on the meteorological site (150 Km is the spatial 
resolution of the 6.6 GHz channel). Air and ground temperature data for the Fall of 1984 were 
obtained from NOAA's National Climatic Data Center in Asheville, North Carolina. Air 
temperature measurements were available for noon and midnight at the meteorologic sites (i.e., 
simultaneously with the satellite pass), but ground temperature measurements were for 7:00 a.m. 
and 7:00 p.m. EST, and were not co-located with the meteorologic sites. Ground temperatures are 
measured at 5 cm depths. Diurnal heating will affect 5 cm temperatures so that there will be 
differences between those temperatures and the effective soil temperatures at the times of satellite 
passes. 


Local area averages at the meteorologic sites were used to define the preliminary boundaries 
in our Freeze Indicator discriminant. For example. Figure 1 illustrates the correlation between 37 
GHz radiobrightness and reported air temperature. The nominal line in these figures is a single 
best fit linear regression in the least squares sense of all local area averages. Individual linear fits 
will differ slightly as shown in Figure 1 (a). We used the nominal line in our discriminant for 
simplicity, but a more sophisticated discriminant might use the actual least squares fit for the local 
area and for the time of day. The discriminant boundaries in Figures 1 (b) and 1 (c) are merely 
estimates based upon the nominal regression and a compromise between midnight and noon air 
temperatures that would imply frozen soil (the lower boundary) and thawed soil (the upper 
boundary). Remember that diurnal temperature gradients will generally cause midnight, sub- 
surface soil temperatures to be warmer than air temperatures, and noon, sub-surface soil 
temperatures to be colder. 

Similarly, local area averages of spectral gradient versus air temperature were the bases for 
the spectral gradient decision boundaries shown in Figure 3. Note that the midnight freeze 
boundary in this example is relatively unambiguous, while a more effective noon freeze boundary 
would be shifted upwards by 0.2 K/GHz. Again, for simplicity in this initial study, we used 
discriminant boundaries that were time and location independent. 

Our 2-parameter Freeze Indicator incorporates the single-band, 37 GHz radiobrightness, 
and a spectral gradient based upon linear regression of 10.7, 18, and 37 GHz radiobrightnesses for 
each 100x100 Km pixel. Based upon the decision boundaries in Figure 1 (b) and 1 (c), the 
likelihood of frozen ground in a 37 GHz pixel, P 37 , is estimated as 


P37 


Tbmax - Tb 37 

Tbmax ■ Tbjnjn 


0 : Tb 3 7 > Tb max 

Tbmin < Tl>37 < Tb 

max 

A u mm 

1 : Tb37 < Tb mm 

\ 

where Tb 37 is the measured 37 GHz radiobrightness, and the preliminary decision boundaries are 


Tbmax = 259 K 
Tbmin = 247 K 


Page 7 


The likelihood of frozen ground based upon spectral gradient decision boundaries in 
Figure 3(a) and 3(b) is p S g, and is estimated as 


Psg ~ 


/ 


0 


[3Tb 1 

3Tb 

1 df /max 

a/ 

i 3Tb i 1 

(3Tb\ 

\ df /max ' 

13/Ln 


3Tb [3Tb \ 

3/ \ 3/ Lax 

[3Tb \ ^3Tb J3Tb| 

" \ 3/ /min 3/ ' 3/ L 


3Tb [ 3Tb 1 
3/ < l 3/ In 


where the preliminary decision boundaries are 
,3Tb, 




= 0.3 K/GHz 


max 


(©) = -0.3 K/GHz 

min 


These boundaries are preliminary in that they were chosen to yield the fewest 
misclassifications in plots of the type shown in Figure 5(a) and 5(b). More refined discriminants 
would incorporate area and time specific decision boundaries. This would be relatively 
straightforward if there were a higher density of weather stations in the test area. As it is, we 
believe that diurnal temperature modeling well yield effective time dependent boundaries, and, 
perhaps, requiring sub-region consistency within a classification will yield improved spatially 
dependent boundaries. The basic sparseness and lack of control of air and ground data should 
prompt some caution about over-interpreting these results. 

Our freeze/thaw discriminant, or Freeze Indicator, is the product of P 37 and psg, and is 
applied at the scale of the 10.7 GHz data. Resolution differences between different frequency 
channels can produce anomalous composite image results if the data were processed at their 
original scale. To avoid these problems, the resolution of the data from each channel is 
compensated to the (coarse) resolution of the lowest frequency channel used in estimating spectral 
gradients (i.e. 10.7 GHz and 100 Km resolution). Under certain constraints upon the 
classification process, these images can be referenced to the higher resolution, 37 GHz format for 
better location of freeze/thaw boundaries (Zuemdorfer, et al., 1989). The effort needed to do this 
would be justified as a part of an improved classification process. Figures 7 through 12 include 
images of the Freeze Indicator for various times during the test period. Black in these images 
indicates a high likelihood of frozen ground. 
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Figure 6(a) and 6(b) show normalized brightness temperatures for midnight and noon, 
respectively, in the northern Great Plains during the Fall of 1984. Normalized brightnesses are the 
average regional brightness at each microwave frequency divided by the average regional air 
temperature. Normalized brightness thus has the dimension of emissivity. Note that there is little 
systematic ordering among the 10.7, 18, and 37 GHz normalized brightnesses during August 
through most of November. However, during the latter half of November through December, the 
normalized brightnesses at midnight are uniformly ordered, 10.7 GHz brightnesses are high, 1 8 
GHz brightnesses are middle, and 37 GHz brightnesses are low. That is, they exhibit negative 
average spectral gradients. The noon normalized brightnesses for December exhibit a similar 
trend, but with exceptions. These are, we believe, illustrations of the law of darkening for frozen 
soils. Soils at midnight in December for the northern Great Plains are very likely to be frozen. 

Performance of the freeze/thaw discriminant is demonstrated in Figures 7-9 where Freeze 
Indicator (FT) images are compared with ground and air temperature measurements for midnight on 
9/20/84, 10/24/84, and 12/9/84. Midnight FI images are shown as better examples of the potential 
of a freeze discriminant. Noon FI images are generally less consistent with meteorologic reports 
because of the contribution of the noontime positive diurnal spectral gradient to the negative frozen 
ground spectral gradient that we discussed in the last section. Areas not covered by the satellite in 
a particular pass are shown in white. Tables 1-3 are summaries of the meteorologic reports. 

On the night of September 20 (Fig. 7), air temperatures throughout the region were near 
60° F and had been above freezing for several days. The FI image shows weak, probably false 
indications of freezing in the prairies of ND, southern Canada, and the rolling glacial terrain east of 
the Red River Valley in Minnesota. While the dry air of the northern prairies permits nighttime 
radiation cooling of the ground to temperatures below that of the air, the more likely explanation 
for the weak freeze indication is short wavelength scattering by the tall prairie grasses in the 
northern great plains, and by woodland areas in Minnesota. However, there are no strong 
indications of freezing in the FI image. 

On the night of October 24 (Fig. 8), air temperatures hovered about freezing throughout the 
area, but had been below freezing at Williston for several days, and generally above freezing 
toward the east (see the temperatures for Fargo, Aberdeen, and Huron in Table 2). The FI image 
shows a strong freeze indication in northwestern ND which is consistent with the temperature 
patterns. Similarly, the definite thaw indication along the Red River Valley is consistent with the 
warmer temperatures reported and the generally more moist soil in the Valley. 

On the night of December 9 (Fig. 9), air temperatures were generally below freezing except 
at Rapid City, SD, and had been below freezing for several days. There was no more than trace 
snow on the ground anywhere in the region. The FI image shows strong freeze indications 
throughout most of the region with weaker indications near Rapid City, and in the Aberdeen-Fargo 
sub-region (Aberdeen is not shown on the December 9 map because its temperature report was 
missing for that date). Again, the FI image is consistent with the temperature record. 

37 GHz radiobrightness and FI image sequences were produced at midnight and noon for 
six-day periods in September, October, and December (Figures 10-12). SMMR coverage is based 
on a 48 hour cycle— midnight (0000 local hours on the date shown), noon (1200 hours on the same 
date), and then midnight again 36 hours later. However, orbit precession causes gaps in the cycle 
and variations in the coverage footprint. Within these constraints, our objective was to observe, if 
possible, weather dynamics reflected in the FI images. 
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The 37 GHz sequence beginning on September 16 (Figure 10 and Table 1) shows the 
moist area associated with the Missouri River, Sakakawea and Devils Lakes in ND, and the 
Missouri River and Lake Oahe in SD. Rain during the night of September 21 appears as a regional 
darkening of the 37 GHz image for midnight on the 22nd. Note that the rain is not picked up in the 
FI image. 

The October sequence (Figure 1 1) is dominated by a cold front passing through the area 
from the northwest with rain and snow beginning on October 19. The region is warmer and drier 
by the 26th. The moisture pattern dominates the 37 GHz image, but only the apparent freeze 
pattern, which generally lags the cold front, is shown in the FI image. Note that strong freeze 
indications follow the cold front but weaken in the south with warming on the 26th. 

The December sequence (Figure 12) is characterized by cold temperatures and snow from 
December 2 through December 5, followed by daytime warming into the 40s (and even 58° at 
Rapid City, SD) by the 9th. The FI images reflect this general coldness, but also show daytime 
thawing toward the end of the period. 


Accomplishments and Open Items 

Freeze Indicator images based upon a preliminary, 2-parameter discriminant— 37 GHz 
radiobrightness and 10.7, 18, and 37 GHz spectral gradient-show relatively good correlation with 
the expected state of moisture in northern Great Plains soils during the Fall of 1984. The 
discriminant is preliminary in the sense that experimental testing of theoretical models needs to be 
done to fully understand the spectral radiobrightness signatures of frozen soils. The concept 
underlying the preliminary discriminant is that frozen soil will exhibit volume scatter darkening at 
shorter microwave wavelengths much like the effect observed in dry snow. Few other phenomena 
cause negative microwave spectral gradients. One such phenomenon is diurnal insolation which 
should cause negative spectral gradients at midnight, but positive spectral gradients at noon. We 
have modeled diurnal insolation and are in the process of tailoring our discriminant to allow for 
diurnal gradients. 

Freeze Indicator images based upon SMMR data effectively map temporal variations in the 
freeze/thaw pattern for the northern Great Plains at the time scale of days. These patterns are 
synchronized with weather patterns, but are not identical. We intend to expand our test data set to 
include several complete seasons. The product would be, in essence, a movie of freeze/thaw 
patterns as weather fronts sweep through the Great Plains throughout several seasons. The 
development of these data from SMMR archives should provide one aspect of hydrologic and 
mesoscale climatic baselines for the region. 

The one significant objective of our initial project that has not been achieved is the 
incorporation of the frozen/thawed soil classification map into regional climate climate models. 
Such classification maps represent new boundary parameters that are not a part of current models. 
A very high priority of the new investigation is to explore how this ffeeze/thaw parameter should 
be incorporated. 

Ancillary to the scientific achievements, but, nevertheless, extremely important to this 
investigation, and to continuing freeze/thaw investigations, are the computational and image 
manipulation tools that have been developed within the Radiation Laboratory partly as a result of 
this project. The primary computational tool is the Fortran coded model for the spectral 
radiobrightness of periodically heated, two-phase media (England, 1989). This model guided our 
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examination of SMMR s spectral response of soil to freezing and thawing, and to diurnal 
insolation. Our image processing tools have evolved from resampling and classification algorithms 
on a general use, Apollo system, to a dedicated image processing system resident in a Sun/4 
Workstation. The Sun/4 system includes two commercial, image processing software packages, 
and peripherals for reading tapes (and, shortly, CD's), digitizing maps and images, and printing 
color images to film and to paper. 



TABLE 1 


Year 

1984 

Sit* 

Air Temp 

CD 

QQ 12 

9/15 

Aberdeen 

40 

63 


Bismark 

39 

64 


Fargo 

46 

63 


Huron 

44 

64 


Miles City 

50 

65 


Rapid City 

44 

63 


Williston 

46 

60 

9/16 

Aberdeen 

46 

65 


Bismark 

53 

67 


Fargo 

47 

69 


Huron 

50 

64 


Miles City 

58 

74 


Rapid City 

52 

70 


Williston 

55 

68 

9/17 

Aberdeen 

57 

71 


Bismark 

56 

83 


Fargo 

57 

73 


Huron 

56 

73 


Miles City 

61 

80 


Rapid City 

64 

84 


Williston 

56 

74 

9/18 

Aberdeen 

64 

81 


Bismark 

57 

80 


Fargo 

64 

81 


Huron 

62 

84 


Miles City 

60 

80 


Rapid City 

63 

90 


Williston 

67 

73 

9/19 

Aberdeen 

59 

64 


Bismark 

55 

88 


Fargo 

65 

90 


Huron 

68 

90 


Miles City 

67 

76 


Rapid City 

60 

90 


Williston 

57 

72 

9/20 

Aberdeen 

62 

65 


Bismark 

58 

63 


Fargo 

54 

62 


Huron 

62 

71 


Miles City 

63 

72 


Rapid City 

59 

67 


Williston 

54 

51 

9/21 

Aberdeen 

59 

82 


Bismark 

58 

81 


Fargo 

59 

77 


Huron 

60 

83 


Miles City 

66 

56 


Rapid City 

60 

83 


Williston 

57 

54 

9/22 

Aberdeen 

61 

65 


Bismark 

52 

60 


Fargo 

72 

66 


Huron 

66 

66 


Miles City 

46 

52 


Rapid City 

45 

60 


Williston 

45 

44 


Cloud 


Snow 

Cover 


Pack 

(x/101 

Prado, (in.) 

(in.) 

QQ 12 

QQ 12 240 Bam 

24h 


2 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

T 

0 

0 

3 

0 

0 

0 

0 

2 

2 

0 

0 

0 

0 

4 

2 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

1 

3 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

10 

0 

0 

0 

0 

0 

3 

- 

- 

0 

0 

3 

6 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

2 

8 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

10 

3 

0 

0 

0 

0 

0 

0 

- 

. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 

- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

• 

• 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

7 

2 
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(a) Midnight . A linear 
regression of the Bismark data 
is shown ("Actual" curve), 
with a regression line for data 
collected throughout North 
Dakota and the surrounding 
region ("Nominal" curve). 


lb) Midnight . Shown with 
the Bismark data is a linear 
regression for data collected 
throughout North Dakota and 
the surrounding region 
("Nominal" curve). Also 
shown are brightness temper- 
ature decision thresholds for a 
frozen or thawed surface. 


(c) Noon . Shown with the 
Bismark data is a linear re- 
gression for data collected 
throughout North Dakota and 
the surrounding region 
("Nominal" curve). Also 
shown are brightness temper- 
ature decision thresholds for a 
frozen or thawed surface. 
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Figl 11 '? 1 37 GHz SMMR brightness temperature versus measured surface air temperature, 

Bismark, North Dakota. Data were collected from 8/1/84 to 12/3 1/84. 
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Figure 2. Soil surface thermal gradient. 
Computed for 10% moist soil at Bismarck, 
ND, for a standard day in September. 
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(a) Midnight 



Ar Temp (K) 


(b) Noon 


Fl fi ur C ?• Frequency gradient versus measured surface air temperature, Bismark, North Dakota 
Data were collected from 8/1/84 to 12/31/84. The frequency gradient is the frequency 
regression slope for simultaneous SMMR brightness temperatures at 37GHz, 18 GHz. 
and 10.7 GHz. Shown with the Bismark data are frequency gradient decision 
thresholds for a frozen or thawed surface. 
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Figure 4 . 


Brightness temperature versus equivalent dry snow. Crater Lake, Oregon, 22 March 
1970 (Edgerton et al„ 1971). Note the short wavelength darkening evident for thick 
snowpacks. 
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Figure 5. Frequency gradient versus SMMR 37 GHz brightness temperature, Bismark, North 
Dakota. Data were collected from 8/1/84 to 12/31/84. Shown with the Bismark data 
are clustering decision thresholds for a frozen, mixed, or thawed surface. Based 
upon ground truth, the solid boxes are frozen, open boxes are thawed, and x's are 
mixed pixels. 
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Figure 6 . 37 GHz, 1 8 GHz, and 10.7 GHz SMMR normalized brightness temperatures versus 

calender day. Measurements were made at irregular intervals from 8/1/84 to 
12/31/84. The normalized brightness temperature of a single SMMR frequency 
channel is the average brightness divided by the average surface air temperature, 
averages are calculated over North Dakota and the surrounding region. 




PAGE 20 


Figure 7 


Air and Ground Temp Night 9/20/84 



AIR TEMP 

Thaw O 
Mixed @ 
Freeze £ 


GNDTEMP 
Thaw | | 
Mixed gg 
Freeze | 



A comparison of reported air and ground temperatures with the 
Freeze Indicator for midnight, September 20, 1984. 
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Air and Ground Temp Night 10/24/84 
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Figure 8. A comparison of reported air and ground temperatures with the 
Freeze Indicator for midnight, October 24, 1984. 
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Air and Ground Temp Night 12/9/84 
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Figure 9. A comparison of reported air and ground temperatures with the 
Freeze Indicator for midnight, December 9, 1984. 
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